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Abstract 

We present the results of the near-infrared (IR) spectroscopy of nine IRAS galaxies (NGC 1266, NGC 
1320, NGC 2633, NGC 2903, NGC 3034, Mrk 33, NGC 7331, NGC 7625, NGC 7714) with the ISLE 
imager and spectrograph mounted on the Okayama Astrophysical Observatory 1.88 m telescope. [Fe II] 
1.257 /im and Pa/3 emission lines were observed for the whole sample while II2 2.121 /im and Br7 lines 
were additionally obtained for two sources, whose flux ratios are used as a diagnostic tool of dominant 
energy sources of the galaxies. We find that the nucleus of NGC 1266 is most likely a low ionization 
nuclear emission-line region (LINER), while NGC 2633 and NGC 2903 possibly harbor active galactic 
nuclei (AGNs). No AGN or LINER signal is found for other objects. In addition, we find the spectral 
features which is indicative of some unusual phenomena occurring in the galaxies, such as the large [Fe 
II] line widths compared to the local escape velocity in NGC 1266. The present work shows the potential 
ability of the ISLE to shed new light on the nature of infrared galaxies, either through a statistical survey 
of galaxies or an exploration of spectral features found in individual objects. 

Key Vifords: galaxies: active — galaxies: evolution — galaxies: individual (NGC 1266, NGC 2903, 
NGC 3034) — infrared: galaxies 
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1. Introduction 

One of the major achievements of the Infrared 
Astronomical Satellite (IRAS) is the discovery of a sig- 
nificant number of galaxies emitting the bulk of their ra- 
diation at far-infrared (IR) wavelengths. Luminous IR 
galaxies (LIRGs, with IR luminosity Lir > IO^^Lq) are 
as numerous as optical starburst and Seyfert galaxies 
with similar bolometric luminosity in the local universe, 
while the more extreme class of ultra-luminous IR galax- 
ies (ULIRGs; LiR, > lO-^^L©) has similar space density and 
bolometric luminosity to optically-selected quasars (Soifer 
et al. 1987). Many studies have been devoted to the ex- 
ploration of the origin of these IR galaxies, and it is now 
becoming the general agreement that strong interactions 
and mergers of gas-rich galaxies are the trigger for the 
majority of them. In the complete sample of IRAS 1 Jy 
ULIRGs, Veilleux et al. (2002) find that nearly 100 % of 
IR galaxies show strong signs of tidal interaction. 

On the other hand, there is still much debate about 
energy sources of IR galaxies. Their IR luminosity can 
certainly come from dust reprocessing of radiations from 
starburst and/or active galactic nuclei (AGNs) activity. In 
this context, it is worth noting that ULIRGs and AGNs 
most probably have evolutionary connection. It is sug- 
gested that major mergers of gas-rich galaxies first form 
a "cool" ULIRG dominated by dusty starburst, which 
is followed by a "warm" ULIRG phase when an AGN 



turns on and starts to heat surrounding dust. Then the 
AGN evolves into the dominant energy source and blows 
away the surrounding dust cocoon, leading to an optically- 
bright quasar (e.g., Lonsdale et al. 2006). Therefore, ob- 
servational investigation of energy sources in IR galaxies 
at various evolutionary stages is a key for understand- 
ing the whole picture of galaxy evolution including AGN 
phase. However, such investigation has been hampered 
by heavy dust obscuration inside IR galaxies. Although 
numerous efforts have been devoted to the subject at var- 
ious wavelengths (e.g., Veilleux et al. 1995; Tran et al. 
2001; Ptak et al. 2003; Farrah et al. 2007; Imanishi et al. 
2007), the results are still largely controversial: different 
studies give the inconsistent estimates of relative contribu- 
tions from starburst and AGNs for the same population. 

In order to address the issue of energy sources in IR 
galaxies, we aim to investigate the presence of AGNs inde- 
pendently from the previous studies, in a large number of 
IR galaxies including ULIRGs, LIRGs, and less luminous 
populations. In this paper, we present the near-IR spectra 
of nine IRAS galaxies with relatively low IR luminosity, 
obtained by the ISLE imager and spectrograph mounted 
on the Okayama Astrophysical Observatory (OAO) 1.88 
m telescope. Near-IR light has a few advantages over 
those at other wavelengths in probing IR galaxies. First, 
it suffers much less from dust obscuration than opti- 
cal light, since the extinction at the former wavelengths 
are only about a tenth of that at the latter wavelengths 
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{Ak ^ O.l Ay). Second, near-IR observations cost much 
less than X-ray, mid- or far-IR observations, hence a 
statistical number of objects can easily be investigated. 
Third, there are many emission lines suitable for probing 
their radiation sources in near-IR spectral region (e.g., 
Matsuoka et al. 2007; Matsuoka et al. 2008). Veilleux et 
al. (1999) demonstrate the power of near-IR spectroscopy 
by observing 39 ULIRGs and finding that at least 50 
% of optically-classified Seyfert 2 galaxies present hid- 
den broad-line-region (BLR) emissions in near-IR lines. 
It means that nuclear regions of many optical Seyfert 2 
galaxies, obscured completely in optical wave bands, are 
actually optically thin at near-IR wave bands. 

We observed the four near-IR emission lines, [Fe II] 
1.257 iim and Pa/3 for the whole sample of nine galax- 
ies and H2 2.121 ^m and Br7 for two of them, whose flux 
ratios have been proposed as a diagnostic tool of domi- 
nant energy sources of emission-line galaxies (Larkin et 
al. 1998; Rodriguez- Ardila et al. 2005). Note that this 
kind of diagnostic diagrams works only with a statistical 
number of sample, since there are always outlying objects. 
In this sense, we do not intend to present conclusive argu- 
ments with the small number of objects observed in this 
work, but rather aim to demonstrate the capability of the 
ISLE to investigate near-IR emission lines of IR galax- 
ies for future projects. Hereafter, [Fc II] 1.257 fim and 
H2 2.121 fim lines are abbreviated to [Fe II] and II2 for 
simplicity. 

2. Observation and Data Reduction 

The target objects were selected from the Imperial 
IRAS-FSC Rcdshift Catalog (IIFSCz; Wang & Rowan- 
Robinson 2009) based on the IRAS Faint Source Catalog 
(FSC). The observability of the four emission lines, [Fe 
II], Pa /3, H2, and Br 7, in the near-IR atmospheric win- 
dow gives a severe restriction on redshifts, while the addi- 
tional constraints come from positions and near-IR bright- 
ness ( J < 12 mag). We summarize the redshifts, the Two 
Micron AU Sky Survey (2MASS; Skrutskie et al. 2006) 
J-band magnitudes, the 60-/im fluxes (/eo/^m) and the IR 
luminosity (^ir) based on the IRAS measurements, and 
the corresponding IRAS names of the targets in Table 1. 
They are not extremely luminous at IR wavelengths com- 
pared to LIRGs or ULIRGs. 

The observation was carried out with the ISLE, a near- 
IR imager and spectrograph mounted on the OAO 1.88 m 
telescope (Yanagisawa et al. 2006; Yanagisawaet al. 2008). 
The instrument uses a HAWAII Ik x Ik array, which pro- 
vides 4.3 X 4.3 arcmin^ field-of-view with a pixel scale of 
0.25 arcsec. The useful wavelength intervals covered by 
the J, H, and K filters are 1.11 - 1.32 fim, 1.50 - 1.79 
/im, and 2.02 - 2.37 fim, respectively. The slit length and 
orientation are fixed to 4 arcmin and east-west direction. 

The observation journal is given in Table 2. For each 
object we first carried out J-band observation, and then 
decided the priority of taking additional H- or A'-band 
data based on the J-band spectrum. As a result, iJ-band 
spectrum was taken for NGC 1266 and AT-band spectra 



were taken for NGC 2633 and NGC 2903 (more details 
are given in the following section). The sky condition was 
sometimes non-photometric during the nights. We used 
the 2.0-arcscc slit, which provides wavelength resolution 
of i? - 1200, 1800, and 1000 at J, H, and K band, re- 
spectively. The choice of the relatively wide slit was effec- 
tive in reducing amplitudes of fringe patterns which often 
dominate measurement errors at red part of the spectra. 
The total exposure times were broken into individual ex- 
posures of 120 sec. and the objects were offset along the 
spatial direction of the slit between adjacent exposures. 
We observed AOV - AOIII stars for fiux calibration im- 
mediately before or after the target exposures at similar 
airmass. 

The data reduction was performed in a standard man- 
ner. After dark subtraction and flat fielding, the sky emis- 
sion was eliminated by subtracting an adjacent offset im- 
age. Then the residual sky background was estimated 
from counts of nearby pixels in the spatial direction, and 
removed. We extracted the target spectra within a 5.0- 
arcsec (20-pixel) aperture, whose size is close to the full 
width at zero intensity of the point spread function, in all 
but the J-band images of NGC 7331 and NGC 1266 re- 
quiring a larger 7.6-arcsec (30-pixel) aperture due to the 
poor seeing. The apertures were centered on spatial peak 
positions of the spectra. All the galaxies have the appar- 
ent sizes of >10 arcsec, hence the resultant spectra sample 
only their central regions. We also note that the esti- 
mated sky backgrounds could contain the contributions 
from the outer galaxies. Such a contamination would re- 
duce the galaxy contributions in the extracted spectra, 
which might results in a little enhancement of the relative 
contribution of AGNs. However, the signal-to-noisc ratios 
of outer galaxy regions in our data are not high enough to 
quantify these contributions. Wavelength calibration was 
achieved by referring to the Ar arc spectra obtained with 
the same instrument configuration as used for the target 
observation. The mean RMS of the calibration are 0.9, 
3.4. and 1.7 A in J, H, and K band, respectively. The at- 
mospheric and instrumental transmission was estimated 
from the observed spectra of the AO standard stars, for 
which we assumed an intrinsic black-body spectrum with 
the effective temperature of 9600 K (Pickles 1998) after 
stellar H I recombination lines were manually removed. 

3. Results and Discussion 

We show the reduced spectra around the four emission 
lines, as well as [Fe II] 1.644 /xm for NGC 1266, in Figure 
1. To the detected lines, we fit Gaussian functions with 
underlying continua represented by tilted lines, i.e.. 



F(v) = ao exp 



(t'-ai)' 



2ai 



+ (as + a^v), 



(1) 



where v represents velocity shift relative to the line cen- 
ters. The free parameters Oq, ai, 02, 03, and 04 are de- 
termined simultaneously by fitting the function to the ob- 
served spectra with the least-x^ method within velocity 
range from -2000 to -^2000 km s'^. We show the fitted 
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Table 2. Observation Journal 
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functions in Figure 1. The equivalent widths (EWs) and 
the full widths at half maximum (FWHMs) of the emission 
lines derived from the best-fit parameters are summarized 
in Table 3. The FWHMs were corrected for the instru- 
mental resolution assuming R = 1200, 1800, and 1000 in 
J, H, and K band, respectively. The upper hmit of 300 
km s^^ was given to the lines with the measured FWHMs 
less than this value. The reduced values are close to 
one in all the cases. Note that we use EWs and flux ratios 
rather than absolute flux values, since the latter is subject 
to the unknown amounts of aperture losses. 

The measured line flux ratios of [Fc II] /Pa (3 and H2/Br 
7 are plotted in Figure 2. Larkin et al. (1998) suggest this 
diagram as a diagnostic tool for energy sources of the line 
emissions. Later the classification scheme is updated by 
Rodriguez- Ardila et al. (2004) and Rodriguez- Ardila et al. 
(2005), who find that AGNs are characterized by the two 
ratios between 0.6 and 2, while the smaller or larger values 
indicate starburst/H II galaxies or low ionization nuclear 
emission-line regions (LINERs), respectively (dotted lines 
of Figure 2). However, note that this diagnostic diagram 
works only with a statistical number of sample, since there 
arc the objects that do not meet the above criteria in the 
compilation of Rodriguez- Ardila et al. (2005). 

The most noticeable object in Figure 2 is NGC 1266. 
Its [Fe II] line flux is at least 10 times larger than Pa /3, 



clearly indicating it is a LINER. Therefore we obtained 
the additional i?-band spectrum and found that another 
[Fe II] line, [Fe II] 1.644 /xm, is also strong in this object. 
The FWHMs of the two [Fe II] lines, -500 km s-\ are 
significantly larger than the local escape velocity (< 340 
km s~^; Alatalo et al. 2011), which is indicative of an en- 
ergetic phenomenon occurring in the galaxy. Actually, the 
presence of an AGN or LINER in the galaxy is implied by 
the radio, optical, and X-ray observations (Alatalo et al. 
2011 and references therein). A powerful molecular wind 
from the nucleus of NGC 1266 is found, which is sug- 
gested to be driven by an AGN since the estimated star- 
formation rate is insufficient for the drive. Our near-IR 
exploration present a new evidence for the LINER nature 
of this galaxy. 

NGC 2633 and NGC 2903 have the consistent [Fe II] /Pa 
/3 and H2/Br 7 ratios with AGNs. H- and if -band spectra 
of NGC 2633 were obtained by Vanzi et al. (1998), who 
find that H2 and Br 7 arc clearly detected while [Fe II] 
1.644 jim is not. The EWs of H2 and Br 7 measured by 
them, 3.0 ± 0.5 A and 5.0 ± 0.5 A, are in excellent agree- 
ment with our results. While Vanzi et al. (1998) classify 
this galaxy as a starburst based on the optical [O I] 0.63 
/.im/H a and near-IR [Fe II] 1.644 /itm/Br 7 line ratios, 
our new data indicate that NGC 2633 may be more like 
an AGN. In reality NGC 2633 sits on the borderline of 
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Fig. 1. ISLE spectra of the IRAS galaxies around [Fe II] 1.257 /im, Pa /3 (observed in J band), H2 2.1213 fim, and Br 7 (in K 
band) emission lines, as well as [Fe II] 1.644 fim (in H band) for NGC 1266. The wavelengths have been converted to the relative 
velocities to the line centers. The target name and the emission line are indicated at the top-left corner of each panel. The fitted 
Gaussian functions are also shown for the detected lines (red dotted lines). The fluxes are given per velocity in arbitrary scale, while 
those of a same object are shown in the same scale. 
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Table 3. Measured Emission-Line Properties 
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Note — 2-tT upper limits are given for the fluxes of undetected emission lines. 
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Fig. 2. Measured lino flux ratios of [Fe II] /Pa /3 and H2/Br 
7. The circles represent the objects with the four lines avail- 
able, while the shaded area (Icr confldence level) or horizon- 
tal lines with arrows (the up and down arrow shows lower 
and upper limit, respectively) represent those without H2 
and Br 7 measurements. The object names are indicated 
near the corresponding symbols. The dotted lines show 
the starburst (SB)/AGN/LINER demarcation proposed by 
Rodn'guez-Ardila et al. (2005), while the dashed lines show 
the approximate envelope of distribution of the objects com- 
piled by them. 



starbursts and AGNs in our diagram and of starbursts 
and "composite" objects in the Vanzi et al. classification 
diagram, suggesting the contributions from both starburst 
and AGN activity. NGC 2903 is a well-known starburst 
galaxy and has been studied at various wavelengths from 
radio to X-ray (e.g., Popping et al. 2010). While the 
majority of the radiation is thought to come from star- 
formation activity, Perez-Ramfrez et al. (2010) raise the 
possibility of a low-luminosity AGN present in this galaxy 
based on its X-ray property. The ncar-IR line flux ratios 
presented here support this possibility. 

The [Fe II] /Pa ^ ratios of NGC 3034, Mrk 33, NGC 
7625, and NGC 7714 are found to be at most 0.3, while 
no observation is available for the measurement of II2 and 
Br 7. In the sample compiled by Rodriguez- Ardila et al. 
(2005), aU the sources with [Fe II]/Pa ^ < 0.3 have Hs/Br 
7 < 0.3. It may imply that the above four objects are most 
likely categorized into starburst galaxies in Figure 2, while 
the additional observations of H2 and Br 7 lines in these 
objects are needed for further discussion on this point. No 
clear evidence for the presence of AGNs has been obtained 
at other wavelengths, despite numerous observations ded- 
icated to these well-known starburst galaxies. Another 
worth-noting point is the P Cygni profile observed in [Fe 
II] and Pa /3 lines of NGC 3034 (see Figure 1). The ab- 
sorption features are detected at '^250 km s~^ blucshift 
relative to the line centers, pointing to the presence of 
outflowing material in front of the line emission region. 
These features are not found in the optical spectra of nu- 
clear star clusters of this galaxy (e.g., Westmoquette et 
al. 2007; Konstantopoulos et al. 2009), which may be due 
to the heavy dust obscuration. Observations with higher 
spatial resolution than the present one are required for 
locating this velocity component in the nuclear region of 
the galaxy. 
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4. Conclusion 

In this paper wc demonstrate the capability of the ISLE 
imager and spectrograph mounted on the OAO 1.88 m 
telescope to investigate near-IR line emissions of IRAS 
galaxies. The observation targets were selected from the 
IIFSCz based on the IRAS FSC. All of them are the 
nearby galaxies at the redshifts z < 0.01. [Fe II] 1.257 fim 
and Pa 13 emission lines were observed in a sample of nine 
galaxies, while II2 2.121 ^m and Br 7 lines were addition- 
ally observed for two of them. Based on the measured line 
flux ratios, NGC 1266 is found to have a LINER, while 
NGC 2633 and NGC 2903 possibly harbor AGNs. On the 
other hand, no AGN signal is found for NGC 3034, Mrk 
33, NGC 7625, and NGC 7714. In addition, we find the 
relatively large [Fe II] line widths of ^^500 km s~^ in NGC 
1266 and the P Cygni profile of [Fe II] and Pa f3 lines in 
NGC 3034. The present work shows the potential ability 
of the ISLE near-IR spectroscopy to shed new light on 
the nature of IR galaxies, either through a statistical sur- 
vey of galaxies on diagnostic diagrams such as Figure 2, 
or an exploration of spectral features found in individual 
galaxies. 
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